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Abstract: In the forthcoming months, the KamLAND experiment will probe the 
parameter space of the solar large mixing angle (LMA) MSW solution as the origin 
of the solar neutrino deficit with z/g's from distant nuclear reactors. If however the 
solution realized in nature is such that Am^^; > 2 ■ 10~^ eV^ (thereafter named 
the HLMA region), KamLAND will only observe a rate suppression but no spectral 
distortion and hence it will not have the optimal sensitivity to measure the mixing 
parameters. In this case, we propose a new medium baseline reactor experiment 
located at Heilbronn (Germany) to pin down the precise value of the solar mixing 
parameters. In this paper, we present the Heilbronn detector site, we calculate the 
Pe interaction rate and the positron spectrum expected from the surrounding nuclear 
power plants. We also discuss the sensitivity of such an experiment to \Ue2,\ in both 
normal and inverted neutrino mass hierarchy scenarios. We then outline the detector 
design, estimate background signals induced by natural radioactivity as well as by in- 
situ cosmic ray muon interaction, and discuss a strategy to detect the anti-neutrino 
signal 'free of background'. 
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1. Introduction 

Over the last years, measurements of atmospheric and solar neutrinos gave increasing 
evidence for massive neutrinos and lepton mixing. Recent results of the SNO exper- 
iment combined with those of Super-Kamiokande showed for the first time directly 
that neutrino flavour conversion occurs [|l], 

Whilst atmospheric neutrino measurements (dominant channel: z/^ — > confine 
the oscillation parameters to 1.4 x 10"^ < ^ml^^ < 4.2 x 10~^ eV^ and sin^ 2Qatm > 
0.9 (90% CL) 0, solar neutrino data (dominant channel: —>■ z/^,t) allow various 
disjoint areas in oscillation parameter space ranging from 10^^'^ to 10^^ eV^ 0, ^, 
They are commonly referred to as small mixing angle (SMA), large mixing (LMA), 
LOW and (quasi-) vacuum (VAC) solution. The main objective of ongoing and 
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upcoming solar and reactor neutrino experiments is to unambiguously identify the 
solution and determine with high accuracy Am^^; and the respective mixing angle. 

Solutions with Am^^; < 10~^ eV^ or small mixing angles can be probed best 
with solar neutrino experiments at sub-MeV energies [0]. Values of Am^^; > 10~^ 
eV^ and large mixing angles can be studied both with long-baseline nuclear reactor 
neutrino experiments |]^, ^, and with sub-MeV solar neutrinos IjTT], [T^, |T3|. Large 



mixing, in particular at large Am^^^ values, is favoured when combining all solar 
neutrino data. The LMA solution gives the best values in global analysis for 
Am^^; = (3.7 - 6.3) x lO'^ eV^ and tan^ ~ 0.35 - 0.38 (sin^ 2Qsoi ^ 0.77 - 0.80) 
[§, ^, l^. Large part of the 10~^ eV^ range is allowed by the current data. It is 



limited by the CHOOZ reactor experiment at ~ 1 x 10 eV |15|. In the following we 



shall focus on the potential of long-baseline reactor neutrino experiments to observe 
neutrino oscillations and to determine Am^^; in case that the LMA solution is realized 
in nature. 

Electron anti-neutrinos (z/g) from nuclear reactors have a continous energy spec- 
trum up to about 10 MeV. They can be detected via inverse beta decay on protons 
VeP for Ey^ > Ethr ~ 1.8 MeV, and their energy is derived from the measured 

positron kinetic energy as Ep^ ~ £'e+ + Ethr- The inverse beta decay cross section, 
including recoil, weak magnetism, and radiative corrections, can be found in [l^ . 



If the distance between nuclear reactor and detector is larger than, or equal to the 
oscillation length, neutrino oscillations become observable as an integral reduction of 
the interaction rate, as well as a periodic modulation of the continous Uf, spectrum. 
For two-neutrino mixing the survival probability is 

with the mixing angle Qsou the mass difference Am^^;, the reactor-detector distance 
L and the neutrino energy E^,^. Since the mixing angle is expected to be large, the 
periodic modulation of the energy spectrum should become clearly visible and the 
value of Am^„; could be derived with high accuracy. 

If however, for a given value of Am^^; the baseline is choosen too long, adjacent 
peaks cannot be resolved experimentally. The shape of the positron spectrum then 
appears unchanged, whilst its normalization is determined by the mixing angle only, 
independent of the actual value of Am^^^. In this case only lower limits on Am^^; 
could be derived. An additional smearing can arise, if several reactors with different 
baseline distances contribute to the interaction rate. 

Currently, there are two experiments which have the sensitivity to explore the 
parameter space of the solar LMA solution with z/e's from nuclear reactors: the 
neutrino signal in the KamLAND experiment in Japan is dominated by nuclear 
power reactors at a distance of 160 km. It is thus sensitive to probe distortions of 
the positron spectrum for values of Am^^, ~ 2 x 10^^ to ~ 2 x 10^'* eV^ [1^, |18|, |l£ 
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The BOREXINO experiment in Italy has a characteristic reactor basehne distance 
of about 750 km, thus would observe spectral distortions for Am^^; ~ 4 x 10~^ to 
~ 4 X 10~^ eV^. Both upper limits are approximate only and depend on the energy 
resolution that will be achieved in the experiments. A resolution of < 10% (a) at 
1 MeV is expected in KamLAND and for BOREXINO an even better resolution 
is expected as the design goal for the light yield is about 400 pe/MeV [|ri|]. 

Several authors pointed out that a dedicated reactor neutrino oscillation experi- 
ment is needed, if Amj.^i > 2 x 10~^ eV^. The authors of Ref. concluded that 
a baseline shorter than that of KamLAND is required in order to determine Am^^; in 
the above quoted range. The authour of Ref. discussed a reactor experiment with 
baseline of a few 10 km to investigate the high Am^ range of the LMA parameter 
space. First experimental details of this study were presented in |^. The authors of 
Ref. ||2T[ studied the achievable accuracy of Am^^^ and tan^B in a hypothetical reac- 
tor experiment with a baseline of 20 km and 3000 events/year. Finally, the authors 
of Ref. p2[ discussed the positron energy spectra for a generic reactor experiment 
with 20 km baseline and the accuracy of Am^^; determination. Moreover, the authors 
pointed out that for a distinct combination of mixing parameters, one can distinguish 
normal from inverted mass hirarchy. In summary, a new reactor neutrino experiment 
would become mandatory in the case that KamLAND observes a reduction of the 
event rate but could not resolve the characteristic oscillation pattern of the positron 
energy spectrum. 

In this paper we propose a location and a design for a reactor neutrino experiment 
which would be dedicated to study, with high sensitivity, neutrino oscillations in the 
HLMA parameter range, i.e. 2 x 10^'* < Am^^^ < 10^^ eV^. Moreover, we 
discuss a new approach to constrain \Ue3\ from the observation of "atmospheric" 
driven oscillation imprinted on the positron spectrum. Finally, we estimate for the 
proposed location and detector design, the background signals induced by natural 
radioactivity as well as by in-situ cosmic ray muon interaction, and discuss a strategy 
to detect the neutrino signal 'free of background'. 



2. Experimental Site 

Requiring the first oscillation peak to appear at about 1 MeV above threshold (i.e 
Ep^ ~ 2.8 MeV), the minimal distance between reactor and detector for Am^^^ = 
2 X 10~^ eV^ given by Eq. (|1.1|) should be approximately 17 km. The detector 
should be located sufficiently deep underground to reduce the cosmic ray muon flux. 
Furthermore, the i/g-interaction rate should be dominated by a single baseline only 
to minimize the smearing of the oscillation pattern from several sources. 

Various location have been investigated over the last year. Our most favoured 
site is a saltmine close to Heilbronn, in the south-west of Germany. More than 2000 
cavernes with approximate dimensions of 15 m x 15 m x 150 m have been excavated 
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at a depth between 180 and 240 m (480 to 640 meter of water equivalent (mwe)). 
Two reactors are in the close, nearly equidistant vicinity: Neckarwestheim (2 cores, 
6.4 GWth) in southern and Obrigheim (1 core, 1.1 GWth) in northern direction. Since 
the saltmine is extended about six kilometers in north-south direction, it is possible 
to select a location with a common baseline of exact 19.5 km to each of the two 
reactors. It is noteworthy that the Obrigheim reactor might be shut down in the 
very near future. In this case the baseline could be selected between 14 and 20 km, 
depending to which parameters the experiment needs to be optimized to. In this 
paper we take as generic location the position with equidistant baselines of 19.5 km 
to the reactors Neckarwestheim and Obrigheim. 



3. Anti-neutrino interaction rate 

The Ue spectrum above detection threshold is the result of /5~ decays of 235,238-[j ^^^^ 
239,241 fission products. Mesurements for '^^^U and 239,24ip^ theoretical calcu- 
lations for 238-[j g^j^g used to evaluate the Ue spectrum [23, 24 ; its overall normalisation 



is known to about 1.4% ||25| and its shape to about 2% [^. As a nuclear reactor oper- 
ates, the fission element proportions evolve in time; as an approximation we use an av- 
eraged fuel composition typical during a reactor cycle corresponding to ^■^^U (55.6 %), 
239pu (32.6 %), (7.1 %) and ^^ipu (4.7 %). The mean energy release (W) per fis- 
sion is then 203.87 MeV and the energy weighted cross section for Uf.p ne^ amounts 
to < cr >fission= 5.825 X 10""^^ cm^ per fission. The reactor power {Pth) is related to the 
number of fissions per second {Nf) by iV/ = 6.241 x lO^^sec"^- (Pt,j[MW])/(iy [MeV]). 
The event rate {Rl) at a distance L from the source, assuming no-oscillation, is then 
Rl = Nf ■ < a >fission -np ■ l/(47rL^), where Up is the number of protons of the 
target. For the purpose of simple scaling, a reactor with a power of 1 GWth induces 
a rate of 447.8 events per year in a detector with 10'^^ protons at a distance of 10 km. 
All relevant European nuclear plants are added in turn in order to compute the Ue 
interaction rate at the Heilbronn site. Five of them (see Tab. P contribute to ~ 92% 
of the total rate Rq; the other reactors contribute less than 1% to the total flux each, 
and 8 % in total. Assuming all reactors running at their nominal power, one expects 
~ 1150 Pe interactions per year in a target containing 10^^ protons (as an example, 
a PXE-based scintillator with a mass of 194 tons contains 10^^ protons). 



4. Neutrino oscillation signatures 

A three-neutrino mixing scenario is required to explain both atmospheric and solar 
anomalies. In that case, flavour eigenstates a = e, h,t and mass eigenstates z = 1, 2, 3 
are related through the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix 
U via the relation z/q, = Yl'i=i ^ai [0; |^ • Assuming a "normal" mass hierarchy 
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Reactor 


Distance [km] 


Power [GWth] 


Rl 


Rl/ Rq 


Neckarwestheim 


19.5 


6.388 


754 


66% 


Obrigheim 


19.5 


1.057 


125 


11% 


Philipsburg 


54 


6.842 


107 


9% 


Biblis 


80 


7.420 


52 


4% 


Grundremmingen 


117 


7.986 


26 


2% 


Others (Europe) 


> 100 


~ 293 


86 


8% 



Table 1: Nuclear reactor data and interaction rate at the Heilbronn site. The 

neutrino interaction rates Rl,Ro are given for a target containing 10'^^ protons with the 
nuclear power plant running at 100% of their nominal power (typical mean values of the 
power vary between 80% and 90%). A PXE-based scintillator (CigHig, molar weight=210 
g/mol) with a mass of 194 tons contains 10^^ protons. 



scenario, mi < m2 < ms, the Ue survival probability can be written |29, 22 



P..^.. = 1 - 2\UMl - \Ue3\') (l - COS (4.1) 
1(1 _ |t/^3|2)2 sin2 20^^ ("l _ COS ^""^1^ 



2' ' ' V 

2|[/e3p(l - It/esH sin' ©12 cos - cos 



\ 2E 2E J 2E 

The first two terms of Eq. ^]l] contain respectively the atmospheric driven [Amli = 
Aml^^) and solar driven (Am^i = Aml^i, 612 ~ Qsoi) contributions, while the third 
term, absent from any two-neutrino mixing model, is an interference between solar 
and atmospheric driven oscillations. We only notice here that Ug^ is the PMNS 
matrix element that couples the heaviest neutrino field to the electron field. The 
relative neutrino interaction rate and spectrum compared with the no-oscillation 
case is given, for a single reactor at a distance L from the detector, by R/Rl = 
J ^{E) ■ Pp^^p^{E) ■ dE, where $ is the z/g production energy spectrum weighted 
by the inverse /3— decay cross section. In what follows, the positron energy spectra 
expected from the five power plants listed in Tab. |l] are added in turn to obtain the 
signal at the Heibronn site. If the KamLAND experiment points out any part of the 
HLMA region as the correct solution of the solar neutrino problem, the experimental 
signature at the Heilbronn site will be a reduction of the rate as well as a distortion 
of the positron energy spectrum. 



4.1 Two- neutrino mixing analysis 



If lUesl vanishes the lepton mixing is radically simplified and the relation 4.1 reduces 



then to ITT], with A 



21 



Aml^i and 612 



Qsoh this approximation has been used 
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by previous reactor neutrino experiments O, |30]. If |f/e3p < 0.01 the two-neutrino 



model remains valid with the small correction sin^ 2Bsoi = (1 ~ |f^e3p)^ sin^ 2Bi2. 
The sensitivity of the experiment depends obviously on the background level that 
can be evaluated from the time variation of the reactor signal; at this stage, assuming 
no background and nuclear power plants running at their full capacity, a detector 
containing lO'^^ target protons would be sensitive to a 10% suppression in rate at a 
3 a level after one year of data taking. Fig. |I] displays the contour lines of equal 
rate suppression, and thus the potential sensitivity at the Heilbronn site over the 
LMA region satisfying Am^^j ^ 3 x 10~^ eV^. In what follows, the finite energy 
resolution of the detector has not been included. The expected positron energy 
spectrum is displayed in Fig. | for several combinations Am^^^ - sin^ 

20soZ- Since 

VfJs with different energies arrive at the Heilbronn detector with different phases, 
shape distortions of the positron energy spectrum are expected. These deformations 
depend strongly on the Am^^; value in the range 1 — 10 x 10^^ eV^; this constitutes 
the main advantage of a medium baseline reactor experiment over other choices, in 
order to measure accurately high Am^^^ values lying in the HLMA area. The sensitity 
to a distorted spectrum could extend up to a few times 10~^ eV^, depending on the 
achievable energy resolution. This overlaps significantly with the parameter space 
probed by the CHOOZ and Palo- Verde experiments [^, 0. Whereas the precision 
measurement of Am^^j depends mainly on the energy resolution, the accuracy of 
the determination of the solar mixing angle will rely on the statistics as well as the 
control of the various backgrounds. 

4.2 Three- neutrino mixing analysis 

Thanks to its simplicity the two-neutrino oscillation model is a powerful tool to 
analyse reactor neutrino data. Nevertheless, on the baseline of interest (~ 20 km), 
assuming Am^^; < Am^^^ and U^i 7^ 0, both i/g oscillations due to ^rnl^i (solar 
driven) and Am^^^ (atmospheric driven) can develop without being averaged. This 
would allow to constrain the relevant |t/e3| parameter which plays a very important 
role in the three-neutrino oscillation phenomenology. For instance, it drives the ^ 
[pji ^ ^e) oscillations of the atmospheric and (very) long baseline experiments. 
At present, the most stringent limit comes from the negative result from the CHOOZ 
experiment^; in the case of interest Am^^; > 2 ■ 10""^ eV^, a three-neutrino mixing 
analysis has constrained |f/e3p < 0.036 at 95 % CL, for Am^^^ = 2.5 x 10~^ eV^ 



and sin Qsoi = 0.27 ||29|. In the forthcoming years, the long baseline accelerator 
experiments aim to reach a sensitivity for |f/e3p down to 0.0015 (see for example 
0); a somewhat weaker constraint on |?7e3p = 0.003 is expected from the K2Rdet 



reactor neutrino experiment [0. On a much longer timescale, a sensitivity improved 



by several order of magnitude seems to be achievable with neutrino factories 



"'^A slighly less stringent constraint has been obtained by the Palo- Verde experiment 
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Improving our knowledge of |?7e3| beforehand is however relevant for the design of 
neutrino factories. For instance, no CP-violation effect would be observable in the 
lepton sector if \Uez\ vanishes or is very small. We note here that all these CP- 
violation effects are unambiguous in long baseline experiment only if Am^^; is known 
with a very high accuracy 



In that context, an additional constraint on |?7e3| could be a by-product of the 
Heilbronn experiment. So as to illustrate the potential sensitivity to \Uez\-, we com- 
pute the positron spectrum obtained at the Heilbronn site with both two and three 
neutrino mixing scenarios, for a few \Uez\ — ^Tn^^i combinations. For simplicity, 
we consider the typical set of LMA parameters Am^^^ = 5 — 80 x 10~^ eV^ and 
sm^2Qsoi = 0.8. Taking into account the SuperKamiokande result on atmospheric 
neutrinos ||^, we fix^ Arn^j^ = 2.5 x 10~^ eV^. The resulting spectrums are re- 
spectively displayed in Fig. |] and Fig. ^ for |f/e3p = 0.04, 0.02. The envelop of the 
positron spectrum is roughly given by the two-neutrino solar mixing, whereas ripples 
are imprinted with the "atmospheric" frequency and an amplitude proportional to 
|f/e3p(l — |f^e3p) (fifst and third term of Eq. |4.1| ). To detect these ripples one should 
be able to resolve two adjacent peaks, which are separated by A{E); this requires 
the energy resolution 6E to satisfy 

with 27r/a = 2.54 ■ ^\m\^^\eV'^\ ■ L\m\. At a baseline of 20 km, for Ej^^ = 3, 4, 5 MeV, 
one has respectively the conditions bE < 0.5, 1.0, 1.7 MeV. This seems to be achiev- 
able according to the energy resolution quoted by CHOOZ (0.40 MeV), and expected 
by KamLAND (FWHM ~ 2.35 x 10%v^ ~ 0.50 MeV) and BOREXINO (FWHM 
~ 2.35 X 7%V^ ~ 0.35 MeV) over the same energy range. In addition, if one does 
not meet the requirement, the constraint on energy resolution could be relaxed by 
~ 25% if one choses the southern site of the Heilbronn mine, located at 14 km of 
Neckarwestheim (which provides the bulk of the z/g flux). 

An adequate statistics is also necessary to detect the "atmospheric" driven rip- 
ples; to obtain a rough estimate of the exposure needed, on can subdivide the positron 
energy spectrum in ~ 0.5 MeV energy bins. The statistical error has to be small 
enough to separate a "bump" from an adjacent "gap" (i.e 2 bins separated by an 
amplitude of about 4 ■ |f/e3p(l — |f^e3p) ), at a /3 sigma level. The expected bin 
content should then satisfy Ub > /lQ\Ue3\^. For jf/egp = 0.04,0.02,0.01 one 
obtains respectively Ub = 40, 160, 640 for /3 = 1 and Ub = 160, 640, 2560 for /3 = 2; 
this leads respectively to ~ 0.9, 3.5, 14 x 10'^^ and ~ 3.5, 14, 56 x 10^^ proton-years 
of exposure (assuming no background, and power plants running full time at their 
nominal capacity). Probing |?7e3p down to 0.01 corresponds to a ~ 2% relative effect 

^Am^(„ is assumed to be known at better than 10%, and slight modifications of it do not change 
notably the results. 
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on average. This is close to the systematic incertainties on the energy spectrum 
production , although those frequencies are expected to be lower than the "atmo- 



spheric" ripple frequencies. Finally, a complete likelihood analysis will be required 
to estimate accurately the detection limit of It/gsl as a function of Am^^;, exposure, 
and backgrounds. 

For completeness, it is worthy of mention that the previous discussion remains 
valid if Am^^; < 2 x 10~^ eV^; indeed, if It/gsl is not too small, whatever the Am^^; 
value measured at KamLAND, an experiment at the Heilbronn site would observe 
a global rate suppression from the "solar" driven neutrino oscillation and "atmo- 
spheric" ripples as discussed previously; this is illustrated in the top-left part of 
Fig. I and Fig. | for Am^^^ = 5 ■ 10"^ eV^ and |f/e3p = 0.04, 0.02. 



As pointed out by Petcov and Piai |22], a medium baseline reactor neutrino ex- 



periment could have the potential to distinguish between normal (mi < m2 < ms) 
and inverted (ms < mi < m2) neutrino mass hierarchy. In this latter case, the 
heaviest neutrino field is mainly coupled to the electron neutrino field. To obtain the 
survival probability from Eq. [4.1| one has to apply the following permutations of in- 



dices: 2^3 and 1 — 2. Eq. ^T!] remains then unchanged (with now Am^g = Am^ 



sol 

and Am^j^ = Am|]^) apart from the amplitude of the interference term that can 
be express as 2|[/eip(l — If/eiP) cos^ ^soh notice that Uf,^ is now constrained if one 
considers the inverted mass hierarchy model. The difference between the two hierar- 
chies, for a medium baseline reactor neutrino experiment, is in principal observable 
if sin^ 26so« 7^ 1, |f^e3/iP > 0.03, and Am^^„ is known with a high precision [|^. At 
the Heilbronn site, the net difference between the two kinds of hierarchy can be seen 
as a shift of the phase of the "atmospheric" driven ripples on the positron energy 
spectrum (Fig. ^). Since the interference term has approximately the same frequency 
as the "atmospheric" driven oscillation term, the energy resolution requirement is 
also given by Eq. |4^ . The exposure needed depends mainly on the |t/e3/i| value as 
discussed previously. The residuals between the two-neutrino and the three-neutrino 
cases normalised to the two- neutrino model (Rate(2i/)-Rate(3z/))/Rate(2z/) are shown 
in the right part of Fig. ^, for a few Am^^^ LMA values and |f/e3/iP = 0.02. For 
Am^^; < 2 ■ 10~^ eV^, the interference term vanishes since Am^^^ — Am^^^ ~ Am^^^ 
and no difference can be seen between normal and inverted hierarchy. Then, Fig. ^ 
indicates that the difference between the two mass hierarchies could be distinguished 
in the narrow range Am^j,; = 2 — 4 ■ 10~^ eV^, while it becomes hardly detectable 
for higher Am^^; values. The ascertainable range is somewhat smaller than the one 
discussed in Ref. with Am^^j = 1 — 5 ■ 10"'* eV^. Indeed, at the Heilbronn site 
the interference pattern is slightly washed out since reactors at different distances 
are involved. 



^The matrix element that couple the electron neutrino to the lightest neutrino field. 



- 8 - 



5. Conceptual detector design 



At this stage, we want to line out the basic concepts of a detector design for the 
Heilbronn site. Detector size and energy resolution depend on the physics goals to 
be addressed: the determination of Am^^^ in the HLMA parameter range can be 
performed with a ~ 100 ton detector, while searches for small oscillation pattern 
driven by the atmospheric Am^^^ with an amplitude proportional to |f/e3p requires 
a high event statistics, and therefore a large target mass. We estimate that ~ 10^^ 
proton-years are required in order to reach a sensitivity of |?7e3p > 0.01. This implies 
a detector similar in size to the KamLAND experiment. In order to resolve adjacent 
oscillation peaks, a light yield of ~ 400 pe/MeV is required for both detector sizes. 
A detailed analysis of the exposure (number of protons x years) vs. sensitivity in 
|?7e3p needs to be carried out carefully. 

The detector concept proposed here is similar to that of the Counting Test Fa- 
cility of BOREXINO ||3^ , however scaled up in size and equipped with a hermetical 
muon veto system. It consists of concentrical spherical volumes with a liquid scintil- 
lator target at its center contained in a transparent vessel, a water buffer surrounding 
it, a layer of photomultiplier tubes (PMT's) viewing the scintillator target, and op- 
tically separated, an outer water buffer equipped with PMT's as a cosmic ray muon 
detector. A schematic view is displayed in Fig. ^. 

Our current best choice for the liquid scintillator is phenyl-xylyl-ethan (PXE, 
CigHig) as a solvent and p-Tp and bis-MSB as fluor and secondary shifter, or alter- 
natively PMP. PXE has a low vapor pressure (< 1.4 x 10~^ hPa at 20 "C) and a flash 
point at 149 °C, and therefore complies with general safety regulations. Moreover, 
the density is 0.99 g/cm^ and thus creates modest boyant forces to the scintillator 
contaiment vessel when emersed in water. Comprehensive studies of PXE based scin- 
tillators were performed in the frame of the BOREXINO project comprising optical 
properties, radioactive impurities and purification methods ^7\, Key re- 



sults include low trace contaminations (< 10"^^g U/g), excellent a—P discrimination, 
and high light yield (310 pe/MeV with 20% optical coverage). 

A Ue event is characterized by a prompt positron event which deposits a visible 
energy between 1 and 8 MeV, followed by a delayed 2.2 MeV gamma event arising 
from neutron capture in hydrogen with r ~ 200 yusec. The minimal energy of 1 MeV 
of the prompt event is due to the positron annihilation in the scintillator. Prompt 
and delayed event are spatially correlated (coincidence volume < 1 m^) and both 
have a /3/7-type pulse shape. This characteristic signature allows to discriminate 
efficiently against backgrounds. 

For the rest of this article, in particular for the estimation of the backgrounds, we 
consider the smallest useful detector size which still allows high precision determina- 
tion of Am^^j in the HLMA parameter range. As generic dimensions we use 300 cm 
radius for the scintillator contaiment vessel, filled with a PXE based scintillator (tar- 
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get mass 112 ton), and an optical coverage of about 30 % providing a photo electron 
yield of 400/MeV. Background estimations, as discussed in the following sections, 
require the PMT's to be located on a sphere with a radius larger than about 500 cm. 
In order to shield against external radiation, an additional external water shielding 
of about 200 cm is necessary. This can be achieved by a cylindrical tank of about 
14 m in diameter and similar in height. A schematic view of the detector layout is 
displayed in Fig. ^. 

As a result of the background considerations we note that a detector can be real- 
ized without using a gadolinium loaded scintillator, since the accidental background 
rate can be suppressed to values well below one event per year. Therefore, we can 
achieve the highest possible light yield with liquid scintillators and hence optimize 
pulse shape discrimination, position reconstruction, as well as energy resolution. 



6. Backgrounds 

Oscillation parameters can only be determined with high accuracy if the backgrounds 
to reactor neutrino interaction are small. In particular, a possible observation of small 
|?7e3| values depends critically on the fraction of background events which survive the 
acceptance cuts. 

Backgrounds from primordial and man-made radioactivity can perturb neutrino 
detection, as well as backgrounds induced by cosmic ray interaction. Previous devel- 
opments and experiments provide quantitative guidance on how to design an exper- 
iment which, in principle, allows to measure Ue's "free of background" . In particular 
the Counting Test Facility of BOREXINO ||, the CHOOZ experiment H and 



measurements of muon induced production of radio-isotopes |^ provide a quantita- 
tive database as well as operational experiences. 

6.1 Non-reactor z/g signals 

Further interferences with reactor neutrinos can arise from geo-physical z/g's. Elec- 
tron anti-neutrinos which are emitted in beta disintegration of the decay products 
of primordial ^ss-y g^^^ 232r|n]^ ^^^q interior of the earth have sufficient energy to 
create background to the reactor neutrino signal. The visible energy of the geo-z/g 
signal is below 2.4 MeV and exhibits a characteristical spectral shape. Following 
Ref. [0, a rate between 5 and 33 events per year and 10^^ protons is expected at the 
Heilbronn site^. The exact value depends on the global abundance of 238]j g^^^ "^^"^Th 
in the earth's crust, mantle and core. First measurements of the geo-Ue flux will be 



performed by the KamLAND and BOREXINO experiments ET], R [TT| and data to 



estimate the interference at the Heilbronn site will be available. Albeit the geo-Ug 



^Heilbronn is located on the continental crust. Therefore, the model developed for the interaction 
rate for the BOREXINO detector can be adopted. 



- 10 - 



rate is expected to be < 3% of the reactor rate, an interference with the reactor 
signal at low energies can not be excluded. A possibility to discard this background 
is to establish an analysis threshold greater than 2.4 MeV. 

6.2 Backgrounds from radioactivity 

Naturally occuring radioactivity can create accidental as well as correlated back- 
grounds. The radon daughters ^^^Bi-^^^Po are a potential background of the latter 
type, because their decay sequence is similar to the v^. tag. However, the fast /5 — a 
decay with r = 237 /isec can be discriminated by the energy deposition of the ^^^Po 
a-decay, quenched to ~ 0.8 MeV (i.e. > lOcr separated from the delayed 2.2 MeV 
7), and in addition by pulse shape information. 

Neutrons from spontaneous fission or (a, ra) -react ions can produce recoil protons, 
subsequently thermalize and capture on hydrogen atoms. Selection of high purity 
materials for detector construction, passive shielding together with pulse shape dis- 
crimination provides an efficient handle against this type of correlated background. 
However, high energy neutrons due to cosmic ray muon interaction may contribute 
significantly to the background. This is discussed in detail in Sec. 6.3.2. 

Gamma and beta signals in the scintillator volume may generate accidentally 
background events which mimick z/g interactions. The accidental background rate 
hacc is given by hacc ~ ^J^dTdYdYdet- Here 6p and are the specific background rates 
(in units of sec~^m~'^) for the prompt and the delayed events, respectively. The 
time window for the coincidence is given by r^, the coincidence volume by V^, and 
Vrfet is the total detetction volume. In order to estimate this background rate in the 
detector we use the rather conservative values of r^^ ~ 1 msec and ~ 1 m^. 

If the accidental background rate should be below hacc ~ 1 and hence negligi- 
ble in comparison to the neutrino interaction rate, the condition for the background 
rates for a Vdet ~ 113 rv? detector then reads h^hd < 3 ■ 10^'^ sec^^m^^. We are 
comparing this limit with possible contributions from different background sources. 



6.2.1 External backgrounds 

As external background we describe contributions due to the detector structure ma- 
terial (i.e. PMT's, light concentrators etc.) and the water shielding outside the liquid 
scintillator. 

The most dangerous source for the external background is the strongly pene- 
trating 2.6 MeV gamma line from ^°^T/ in the Thorium chain. The typical Thorium 
activity reached in selected materials for all PMT's and light concentrators in the 
solar neutrino experiment BOREXINO is about ~ 10^ Bq |Tl]. There a coverage of 
30% at a distance of about 6.5m to the center of the detector is achieved. Assuming 
to use PMT's and light concentrators with the same specific activity at a coverage 
similar to BOREXINO, the background scales with the surface of the PM-sphere. 
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The shielding S of 2.6 MeV gammas due to the water in a spherical geometry 
can be parameterized by S* ~ 10~°, with a ~ (i? — r)/0.5m, where R and r are the 
radii of the PM-sphere and the scintillator vessel, respectively. 

Monte-Carlo calculations show that signals due to the external 2.6 MeV gamma 
line contribute with about 40 % to bp and to about 10 % to 6^. The background 
rate in the total scintillator volume can be estimated to be: bp^tot ~ 10^ sec~^ x 
(i?/6.5m)^ X 5*. For R ~ 5.0m and r ~ 3m one obtains bp-tot = 2.4 x 10^^ sec~^ and 
bd-tot = 6 X 10^'^ sec'^. Due to the self-shielding of the scintillator these events are 
distributed in an outer region (roughly between 2.5m and 3m) with a total volume 
of about 50m'^. Hence the specific activity at this geometry (i.e. radius of the PM- 
sphere at 5m) is bpbd ~ 6 ■ 10~^ sec~^m~^ and therefore well below the limit given 
above. 

External backgrounds due to radioactive impurities in the water shielding have 
been estimated for Uranium, Thorium, and Potassium in Monte-Carlo calculations. 
For the prompt background all events with energies above 1 MeV have been counted. 
The energy window for the delayed 2.2 MeV event was chosen to be 3 sigma at a 
light yield of 400 photo-electrons per MeV energy deposition, corresponding to a 30% 
optical coverage for the PMT's and light concentrators. 

The critical concentration limits obtained for Th and U in the water are in the 
range of about 10~^^ g/g 10~^^ g/g, and for K between 10~^° g/g and 10~^^ g/g, 
respectively. The latter obviously gives no contribution to bd, as the 1.46 MeV 
gammal line from ^^K is far below the 2.2 MeV neutron capture energy. 

For Thorium the counting rates in the total volume at a concentration of 5 ■ 
10~^^ g/g iia water are bp-tot ~ 1.5-10"^ sec~^ and bd-tot ~ 3-10^'^ sec~^, respectively. 
As these events are distributed mainly in an outer scintillator volume of about 20m^, 
the specific background therefore is bpbd ~ 10^'^ sec~^ and below the allowed 
level given above. 

For Uranium the allowed numbers are very similar and one can conclude, that 
for both, U and Th, a concentration lower than 10~^^ g/g sufficient. 

These requirements for the water purity have been achieved in the Counting Test 
Facility (CTF) of the BOREXINO experiment p9|. Special care has to be taken for 
Radon in the water, which is often found to be out of the equilibrium of the U-chain. 
The limit given above translates to a tolerable specific Rn- activity of ~ 10^ mBq/m?. 
However, also this (and even better) values have been achieved in the CTF. 



6.2.2 Internal backgrounds 

The acceptable single counting rates for U, Th and K within the liquid scintillator 
were estimated in Monte-Carlo calculations. Also here a 3 sigma window around the 
2.2 MeV line for the delayed event was chosen. The results show, that the limits for 
concentrations of U and Th are in the level of about 10^^^ q/q. 
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The results in detail: at 1 ■ 10^^^ g/ g for Uranium hp^tot ~ 1 ■ 10^^ sec^^ and 
bd-tot ~ 2 ■ 10~^ sec~^, are the total counting rates in the scintillator volume. As 
these events are distributed homogeneously all over the scintillator volume (i.e. ca. 
lOOm^), the specififc counting rates combine to hpbd ~ 2 -10"^ sec~^ m~^, at the design 
goal. A dominant part of this background comes from ^^^Bi decay. Since the ^^^Bi- 
^^^Po coincidence can be tagged, one can reduce the accidental backgound further 
by a about a factor 10. This translates into a a limit for the trace contamination of 

< 1 • 10-12 g\]/g. 

The same limit is basically valid for Thorium. As again the K-activity does not 
contribute to the delayed signal, the corresponding concentration limits are higher 
by about 2 orders of magnitude, i.e. in the range of about 10"^^ g/ g. 

In the CTF and with neutron activation analysis of BOREXINO it has been 



shown that even better limits have been reached with PXE as liquid solvent and 
p-TP as wavelengthshifter. The measured values of U and Th of the scintillator as 
achieved from the company was already in the 10"^^ g/g range. After purification 
in a Silica-Gel column and using water extraction limits in Th and U of 2 ■ 10^^^ g/g 
and 1 ■ 10"^^ g/g have been obtained, respectively. In Tab. |] we summarize the 
requirements for the principal detector components. 



Source Requirements 

PMT's and concentrators distance > 500 cm to center 

Water shield < 5 x lO'^^ g/g, U,Th 

< 5 X 10-13 g/g, K 
Liquid Scintillator < lO-^^ g/g, U,Th 
< 10-" g/g, K 

Table 2: Summary of background requirements due to radioactive trace impu- 
rities. Trace impurities of the PMT's and light concentrators have been taken as obtained 
in the BOREXINO experiment |11|. The limit stated for the liquid scintillator impurities 



does not include the possibility to discard the ^^^Bi-^i^Po coincidence on an event by event 
basis. 

6.3 Backgrounds induced by cosmic ray muons 

The rock overburden of the Heilbronn site is 180 m at "Kochendorf" , the northern 
area of the Heilbronn saltmine, and 240 m at "NeuboUinger Hof" , the south western 
area. With a standard rock density of 2.65 g/cw? the depths correspond to about 
480 m and 640 m of water equivalent (mwe). The underground muon fluxes, its 
zenith angle distribution and mean energies can be calculated following Ref. [42]. The 
vertical intensity at a depth of 480 (640) mwe is 9.8 (4.7) x 10-^ cm-^ s"i sr"i and the 
total flux, taking into account the flat surface topology, is 2.1 (1.0) x 10~^ cm^^ s-^ or 
760 (360)/h/m2. The mean energy of the vertical flux corresponds to 72 (110) GeV. 
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Cosmic ray muons will be the dominating trigger rate at this depth. About six 
muons per second will cross a spherical scintillation target with a volume of 113 
m^. The energy deposition corresponds to about 2 MeV per cm path length which 
provides a strong discrimination tool. 

6.3.1 Muon induced production of radioactive isotopes 

Long-lived muon induced isotopes can not be correlated to the primary muon in- 
teraction if their life time are much longer with respect to the characteristic time 
between two subsequent muon interactions. They can contribute significantly to the 
single trigger rate as well as simulate the v^. coincidence signature. 

Cross sections of muon induced isotope production on liquid scintillator targets 
(^^C) have been measured by the NA54 experiment at the CERN SPS muon beam 
at 100 GeV and 190 GeV muon energies [|^. The energy dependence was found 
to scale as (Ttot{E^) oc E^" with a = 0.73 ± 0.10 averaged over the various isotopes 
produced. 

Cross sections and interaction rates are calculated for a depth of 480 mwe and 
summarized in Tab. |3|. 





Isotopes 


Tl/2 


Eraax 

(MeV) 


(7(72GeV) 
(/ibarn) 


(sec^^) 






12b 


0.02 s 


13.4 


n.m. 








"Be 


13.80 s 


11.5 


< 0.96 


< 8.7 X 10^5 






"Li 


0.09 s 


20.8 


n.m. 








9Li 


0.18 s 


13.6 


1.0 ±0.4* 


(9.2 ±3.5) X 10^ 


-5 




^Li 


0.84 s 


16.0 


2.3 ±0.9 


(2.1 ±0.8) X 10" 


-4 




8He 


0.12 s 


10.6 


1.0 ±0.4* 


(9.2 ±3.5) X 10" 


-5 




^He 


0.81 s 


3.5 


8.0 ± 1.6 


(6.9 ± 1.7) X 10" 


-4 


(3+, EC 


"C 


20.38 min 


0.96 


453 ± 81 


(4.0 ±0.8) X 10" 


-2 




IOC 


19.30 s 


1.9 


61 ± 10 


(5.5 ±0.9) X 10- 


-3 






0.13 s 


16.0 


2.4 ± 1.2 


(2.1 ± 1.1) X 10" 


-4 




^B 


0.77 s 


13.7 


3.3 ± 1.0 


(2.9 ±0.9) X 10" 


-4 




^Be 


53.3 d 


0.478 (EC, 7) 


100 ± 20 


(9.2 ± 1.7) X 10" 


-3 



Table 3: Radioactive isotopes produced by muons and their secondary shower 
particles in liquid scintillator targets. The rate is given for a target of 5.1 xlO^o ^"^C 
(with a C/H ratio of 16/18 in a pxe-based scintillator this corresponds to a 112 ton target) 
at a depth of 480 mwe. Because of the positron anihilation the visible energy in (5'^ decays 
is shifted by 1.022 MeV. *: ^Li and ^He decays could not be separated experimentally. 
Therefore, the cross section given, corresponds to the sum of ^Li and ^He production. 

The dominating background for muon induced single events comes from the decay 
of ^^C with a rate of about 3.5 x 10^ per day or 3.6 x 10~^s^^ m~'^ PXE scintillator. 
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It thus suffices the constraint on accidental background. Correlated backgrounds 
induced by muon interactions, can be created by /?-neutron instable isotopes. Of 
this type of events, with ^^C as target, only ^He, ^Li and ^^Li are produced. The 
cumulative production of *He and ^Li together amounts to about 3 x 10^ per year for 
our generic target size. The characteristic signature of this class of events consists of a 
four- fold coincidence {fi ^ n ^ (3 ^ n). The initial muon interaction is followed by 
the capture of spallation neutrons within about 1 ms. The time scale of the /3-decay 
of the considered isotopes is on the order of a few 100 ms (c.f. Tab. ^) again followed 
by a neutron capture. This signature allows one to discard with high efficiency this 
class of events maintaining a high acceptance probability for real Ue events. 

6.3.2 Muon induced neutrons 

Muon induced neutron production can be estimated from the results of the CTF 
experiment At LNGS the muon mean energy is 320 GeV and the flux corre- 



sponds to 1.16 m ^h ^. Scaling the CTF production rate of 0.3 d ^t ^ to a depth 
of 480 mwe, we obtain about 7.2 x 10^ d^^. About a factor three higher is estimated 



following the compilation of Ref . . A muon preceeding a neutron, therefore must 
not to be mistaken for a positron event with an efficienecy better than one in 10^ 
events. 

A further source of background are neutrons which are produced in the sur- 
rounding rocks by radioactivity and in cosmic ray muon induced hadronic cascades. 
Neutrons can create a recoil proton and be captured by hydrogen atoms after ther- 
malization mimicking a Ue event. At depths below several tens of meter, radioactive 
decays followed by {a,n) reactions dominate the overall flux. However, the energies 
involved are low and a water buffer of a modest thickness efficiently shields these 
neutrons. High energetic neutrons are produced in nuclear cascades up to GeV ener- 
gies and above. Since the primary cosmic ray muons is not penetrating the detector, 
they are invisible and can not be used for discrimination purpose. High energetic 
neutron fluxes have been estimated for the Gran Sasso Laboratory and amount to 
about 25 m~^ year~^ |^3|]. Scaling to the depth of the Heilbronn site gives a flux 
of 5 X 10^ m~^ year~^. Assuming a surface of 113 m^ an overall reduction factor 
of about 5 X 10^ is required. This challenging requirement needs to be obtained by 
passive shielding together with active recognition of the recoil nuclei by pulse shape 
discrimination techniques. 



7. Conclusions 

In this paper we present a possible new reactor neutrino experiment with a baseline 
of ~ 20 km, dedicated to the accurate determination of the solar mixing parameters 
in the HLMA range Am^^j ^ 2 x 10^^ eV^. In this parameter range the KamLAND 
experiment would observe only a suppression of the total rate, while for the proposed 
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experiment the oscillations would become directly visible as an unique pattern in 
the positron energy spectrum. High precision determination of Am^^; would then 
become feasible. Moreover, this would provide essential input for future long baseline 
neutrino experiments. 

We showed that an experiment at the Heilbronn site with a minimal detector 
mass of about 100 tons could be realized yielding an event rate of about 5 x 10^ 
per year (in the case of no-oscillation) at a very low background rate of less than 
one event per year. Furthermore, we discussed that small effects related to |?7e3| 
could be investigated if one increases the detector mass to about 1 kton. |f/e3p could 
indeed be constrained down to about 0.01 and the normal and inverted neutrino 
mass hierarchies could be separated for a distinct combination of parameters. 
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Figure 1: Contour lines of equal suppression by rates only. The first left curve 
represents a 10% rate suppression of the v'g rate, while each other line marks an additional 
10% suppression step. Assuming no background, and nuclear power plants running at their 
full capacity, one can detect a 10% rate suppression at a 3 a level after one year of data 
taking with a detector containing 10^^ target protons. 
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Figure 2: Positron energy spectra for various oscillation solutions assuming 
the power reactors to operate at their nominal power and a detector mass containing 10'^^ 
protons. The spectrum is binned in 250 keV intervals. The finite energy resolution has not 
been included. The black dotted line corresponds to the case of no oscillation. The thick 
light gray (green) line corresponds sin^ '^.Qsoi = 0.6, the medium thick black line (blue) to 
sin^ 2Qsoi = 0.8 and the thin dark gray (red) line to sin^ 2Qsoi = 1.0. 
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Figure 3: Difference of the positron spectrum simulated at the Heilbronn site, 
considering the 2u and mixing scenarios for various oscillation solutions. 

The spectrums are not binned and the finite energy resolution is not included. The thick 
(black) dotted line corresponds to the no-oscillation case. The thick (black) line corresponds 
to the 2v oscillation case (the small correction due to \Ue-i\ is included in the effective 
mixing angle). The (black) dotted line and the thin (green) dashed line represent the 2>v 
oscillation case, accounting respectively for the normal and inverted hierarchy scenario. 
The parameter choosen are sin^ 29soi = 0.8, Am^^^ = 2.5 x 10^^ eV^, and |C/e3p = 0.04 
(|C/eip for the inverted hierarchy). 
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Figure 4: Three neutrino mixing positron spectrum simulated at the Heil- 
bronn site, for various oscillation solutions (left) and corresponding residuals 
(Rate(2z^)-Rate(3z^)/Rate(2z^) (right). The features of these graphs are the same as 
those of Fig. ^, except that |J/e3p = 0.02 (|J7eiP for the inverted hierarchy). 
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Figure 5: Schematic view of the proposed detector. Depicted is the minimal size 
of a detector to study the HLMA region with high accuracy at the Heilbronn underground 
site. About 112 tons of a PXE based liquid scintillator is contained in a transparent sphere 
surrounded by ultra-pure water as a passive shielding. The design goal is to achieve a 
light yield of about 400 pe/MeV which requires an optical coverage of about 30% provided 
by the surrounding PMT's and light concentrators. The PMT's are mounted on an open 
structure which separates optically the outer part of the detector, used as a muon veto. At 
this site, the depth varies between 480 and 640 mwe reducing the muon flux to about 76 
to 36/h/m^. The Ue interaction rate would be about 5 x 10^ per year in the no-oscillation 
case. The detector design allows the background rate to be less than one event per year, 
provided that the design specifications of the scintillator and water buffer are met. To 
investigate |?7e3p down to ~ 0.01 a neutrino target of about 1 kt is needed which implies to 
enlarge the scintillator sphere to 6.2 m radius. The tank diameter as well as height would 
then scale to about 20 m. 
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